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Abstract. We present results from an SED analysis of two lensed high-z objects, the z = 6.56 
galaxy HCM6A behind the cluster Abell 370 discovered by Hu et al. (2002) and the triple arc at 
z ~ 7 behind Abell 2218 found by Kneib et al. (2004). For HCM 6A we find indications for the 
presence of dust in this galaxy, and we estimate the properties of its stellar populations (SFR, 
age, etc.), and the intrinsic Lya emission. From the "best fit" reddening {E(B — V) ~ 0.25) 
its estimated luminosity is L ~ (1 — 4) x 10 11 Lq, in the range of luminous infrared galaxies. 
For the arc behind Abell 2218 we find a most likely redshift of z ~ 6.0-7.2 taking into account 
both our photometric determination and lensing considerations. SED fits indicate generally a 
low extinction but do not strongly constrain the SF history. Best fits have typical ages of ~ 3 to 
400 Myr. The apparent 4000 A break observed recently by Egami et al. (2004) from combination 
of IRAC/Spitzer and HST observations can also well be reproduced with templates of young 
populations (~ 15 Myr or even younger) and does not necessarily imply old ages. Finally, we 
briefly examine the detectability of dusty lensed high-z galaxies with Herschel and ALMA. 



1. Introduction 

Little is known about the stellar properties, extinction, and the expected intrinsic 
Lya emission of distant, high redshift galaxies. Indeed, although it has in the recent 
past become possible through various techniques to detect already sizeable numbers of 
galaxies at z > 5 (see e.g. the reviews of Taniguchi et al. 2003 and Spinrad 2003) the 
information available on these objects remains generally scant. For example, in many 
cases the galaxies are just detected in two photometric bands and Lya line emission, 
when present, serves to determine the spectroscopic redshift (e.g. Bremer et al. 2004, 
Dickinson et al. 2004, Bunker et al. 2004). Then the photometry is basically used to 
estimate the star formation rate (SFR) assuming standard conversion factors between 
the UV rcstframe light and the SFR, and nothing is known about the extinction, and the 
properties of the stellar population (such as age, detailed star formations history etc.) 

At higher redshift (z > 6) even less information is generally available. Many objects are 
found by Lya emission, but remain weak or sometimes even undetected in the continuum 
(e.g. Rhoads & Malhotra 2001, Kodaira et al. 2003, Cuby et al. 2003, Ajiki et al. 2003, 
Taniguchi et al. 2004) . In these cases the Lya luminosity can be determined and used to 
estimate a SFR using again standard conversion factors. Also the Lya equivalent width is 
estimated, providing some possible clue on the nature of these source. However, this has 
lead to puzzling results e.g. for the sources from the LALA survey (Malhotra & Rhoads 
2001, Rhoads et al. 2003) leaving largely open the question of the nature of these objects, 
their stellar populations, extinction etc. 

Strong gravitation lensing is extremely "helpful" for a large number problems discussed 
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at this conference, including also the present one. In particular strong lcnsing has allowed 
to detect several of the highest redshift galaxies known today (e.g. Ellis et al. 2001, Hu 
et al. 2002, Kncib et al. 2004, Pello et al. 2004a, and the review of Pello et al. 2003). 
Also, thanks to the lcnsing magnification, it has been possible to obtain photometric 
observations of reasonable quality in several bands for some of these objects. For exam- 
ple it has even very recently been possible to image a z ~ 7 galaxy with the Spitzer 
observatory at 3.6 and 4.5 /an (Egami et al. 2004) ! As we'll show below (Sects. 2 and 3) 
this allow us to perform a quantitative SED analysis to constrain properties of the stellar 
populations, such as age and star formation (hereafter SF) history (burst or constant 
SF?), their extinction, intrinsic Lya emission etc. A detailed account of this work will be 
published elsewhere (Schaerer & Pello 2004). 

As such, gravitational lcnsing provides a unique opportunity to learn more about some 
selected high-z galaxies. If generalised and applied to larger samples in the near future, 
systematic studies of the properties of lensed high-z galaxies could provide unique insights 
and complementary information to other dccp/ultra-deep surveys targetting blank fields. 
Also, extensions to wavelengths beyond the optical and near-IR with existing facilities 
(e.g. in the radio, mm, and possibly sub-mm) and future observatories should be of great 
interest, as briefly outlined for Herschel and ALMA in Sect. 4. 



2. Stellar populations and dust in a lensed z = 6.56 starburst galaxy 

The lensed z = 6.56 galaxy HCM6A was found by Hu et al. (2002) from a narrow-band 
survey in the field of the lcnsing cluster Abcll 370. Its redshift is established from the 
broad-band SED including a strong spectral break, and from the observed asymmetry of 
the detected emission line identified as Lya. 

We have recently analysed the SED of this object by means of quantitative SED fitting 
techniques using a modified version of the Hyperz code of Bolzonclla ct al. (2000) f . The 
observed VRIZJHK' data are taken from Hu et al. (2002). The gravitational magnifica- 
tion of the source is fj, = 4.5 according to Hu et al. The main free parameters of the SED 
modeling are the spectral template, extinction, and the reddening law. Empirical and 
theoretical templates including in particular starbursts and QSOs (SB+QSO templates), 
and predictions from synthesis models of Bruzual & Chariot (BC+CWW group) and 
from Schaerer (2003, hereafter S03) are used. 

Overall the SED of HCM 6 A (see Fig. 1) is "reddish", showing an increase of the flux 
from Z to H and even to Kf. From this simple fact it is already clear qualitatively that 
one is driven towards stellar populations with a) "advanced" age and little extinction 
or b) constant or young star formation plus extinction. However, for HCM6A a) can be 
excluded as no Lya emission would be expected in this case. 

Quantitatively, the best solutions obtained for three "spectral template groups" are 
shown in the left panel of Fig. 1. The solutions shown correspond to bursts of ages ~ 
50-130 Myr and little or no extinction. However, as just mentioned, solutions lacking 
young ( < 10 Myr) massive stars can be excluded since Lya emission is observed. The 
best fit empirical SB+QSO template shown corresponds to the spectrum of a metal-poor 
starburst galaxy with an extinction of Ay ~ 1. On the basis of the present observations 
a narrow line (type II) AGN cannot be ruled out. To reconcile the observed SED with 

f To convert the observed/adjusted quantities to absolute values we adopt the following 
cosmological parameters: Q m = 0.3, J7a = 0.7, and Ho = 70 km s _1 Mpc -1 . 

J The significance of a change of the SED slope between JH and HK seems weak, and difficult 
to understand. 




Figure 1. Best fits SEDs to the observations of Abell 370 HCM 6A. The red crosses indicate 
the corresponding model broad band fluxes. Solid lines show the best fit for a template from 
the BC+CWW group, dotted from SB+QSO group, and dashed from the S03+ group Left: 
Observed spectral range. Right: Predicted SED in Spitzer/IRAC domain for best fit models. 
Dashed lines show the bursts from the BCCWW and S03+ template groups. The dotted line lis 
the spectrum of SBS 0335-052 from the SB+QSO group with additional Av = 1. The solid lines 
show best fits for constant star formation using different extinction/attenuation laws (Calzetti 
starburst law versus SMC law). The solid triangles illustrate the IRAC point-source sensitivity 
(lir) for low and medium backgrounds excluding "confusion noise". 



Lya, a young population or constant SF is required. In any of these cases fitting the 
"reddish" SED requires a non negligible amount of reddening. 

Although all best fit models require reddening, this result is at present indicative and 
need to be firmed up. Quantitatively (e.g. for constant star formation, solar metallicity 
models, Calzetti law) Ay is typically ~ 0.5-1.8 mag at the 68 % confidence level. Also 
somewhat smaller extinction can be obtained if the steeper SMC extinction law of Prevot 
et al. (1984) is adopted. Zero extinction cannot be ruled out at the ~ 2a level. Better 
photometric accuracy, especially in the JHK bands, is needed to reduce the present 
uncertainties and hence confirm the indication for dust. 

From the best fit constant SF models we deduce an extinction corrected star formation 
rate of the order of SFR(UV) ~ 11 - 41 M© yr" 1 for a Salpeter IMF from 1 to 100 M 
or a factor 2.55 higher for the often adopted lower mass cut-off of 0.1 M©. For continuous 
SF over timescales tsF longer than ~ 10 Myr, the total (bolometric) luminosity output 
is typically ~ 10 10 L Q per unit SFR (in M© yr" 1 ) for a Salpeter IMF from 1-100 M©, 
quite independently of metallicity. The total luminosity associated with the observed SF 
is therefore L ~ (1 — 4) x 10 11 L©, in the range of luminous infrared galaxies (LIRG). For 
t SF ~ 10 Myr the estimated stellar mass is M* « t SF x SFR - (1 - 4) x 10 8 M©. Other 
properties such as the "Lya transmission" can also be estimated from this approach. A 
relatively high Lya transmission of ~ 20-50 % but possibly up to ~ 90 % is estimated 
from our best fit models (see Schaerer & Pello 2004). 

It is interesting to examine the SEDs predicted by the various models at longer wave- 
lengths, including the rest-frame optical domain, which is potentially observable with the 
sensitive IRAC camera onboard the Spitzer Observatory and other future missions. In 
the right panel of Fig. 1 we plot again the 3 best fits. We see that these solutions have 
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fluxes comparable to or above the detection limit of IRAC/Spitzer f. On the other hand 
the strongly reddened constant SF or young burst solutions do not exhibit a Balmcr 
break and are hence expected to show fluxes just below the IRAC sensitivity at 3.6 /jm 
and significantly lower at longer wavelengths. As Lya emission is expected only for the 
reddened SEDs the latter solutions are predicted to apply to HCM 6A. If possible despite 
the presence of other nearby sources, IRAC/Spitzer observations of HCM 6A down to 
the detection limit or observations with other future satellites could allow to verify our 
prediction and therefore provide an independent (though indirect) confirmation of the 
presence of dust in this high-z galaxy. 



3. A lensed galaxy at z ~ 6 7 behind Abell 2218 

This interesting triply imaged object, a possible z ~ 7 galaxy, has recently been 
discovered by Kneib et al. (2004, hereafter KESR) from deep Z band observations 
with ACS/HST. In the meantime it has also been observed with Spitzer (see Richard 
et at, these proceedings; Egami et al. 2004). The currently available observations in- 
clude V606W (undetected) , Ismw, z 850 lp, J, Hn w, Hi 60 \y, and 3.6 and 4.5 fim with 
IRAC/Spitzer. The photometry from these authors has been adopted here to analyse the 
properties of this object in a similar way as for HCM 6A. In practice, small differences 
are found in the published photometry; we therefore adopt three different SEDs (SED1- 
3) to describe this object (see Schaerer & Pello 2004 for details). No emission line has 
so far been detected for Abell 2218 KESR. Its spectroscopic redshift remains therefore 
presently unknown but the well-constrained mass model for the cluster strongly suggests 
a redshift z ~ 6.5-7 for this source. The magnification factors of both images a and b is 
[i = 25 ± 3, according to KESR. 

As a spectroscopic redshift has not been obtained (yet) for this galaxy we here exam- 
ine its photometric redshift estimate. In Fig. [21 (left) we show the photometric redshift 
probability distributions P(z) for the three SEDs (SED1-3) of Abell 2218 KESR using 
three spectral template groups and adopting a minimum photometric error of 0.15 mag. 
For each redshift, P(z) quantifies the quality of the best fit model obtained varying all 
other parameters (i.e. extinction, /Lyf , spectral template amoung template group). Given 
the excellent HST (WFPC2, ACS and NICMOS) photometry, P{z) is quite well defined: 
the photometric redshift ranges typically between z p hot ~ 5.5 and 7.3. Outside of the 
plotted redshift range P(z) is essentially zero. 

To summarise (but cf. Schaerer & Pello 2004), given the absence of a spectroscopic 
redshift, a fair number of good fits is found to the observations of Abell 2218 KESR when 
considering all the free parameters. Three of them are illustrated in Fig. [3 (right). The 
main conclusions from these "best fits" are: 

1 ) Generally the determined extinction is negligible or zero quite independently of 
the adopted extinction law. For few empirical templates we find good fits requiring an 
additional Ay ~ 0.2-0.6 mag, depending on the adopted extinction law. 

2) Although generally burst models fit somewhat better than those with constant star 
formation among the theoretical templates, the data does not strongly constrain the star 
formation history. 

3) Typical ages between ~ 15 and 400 Myr are obtained. A reasonable 1-a upper bound 
on the age of ~ 650 Myr can be obtained assuming constant star formation. Young 
solutions (~ 15 and even younger) are obtained with burst models or some empirical 



f See http : //ssc . spitzer . caltech . edu/irac/sens . html 
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Figure 2. Left: Photometric redshift probability distributions P(z) of Abell 2218 KESR using 
three spectral template groups (solid, dotted, long dashed), three different variants of the SED, 
and the photometry from all filters in which the object is detected (Igi4w to 4.5 iim) assuming 
a minimum photometric error of 0.15 mag. Right: Best fits SEDs to the observations of Abell 
2218 from Egami et al. 2004 (SED2) The red crosses indicate the corresponding model broad 
band fluxes. The solid line shows the best fit for a template from the S03+ group, and dotted 
from the SB+QSO group. The redshift for these solutions are z ~ 6.63 and 6.54 respectively. 
See text for more information 



templates. The relatively modest strength Balmer break observed between the HST and 
Spitzer broad-band photometry does not necessarily imply old ages. 

4) Given degeneracies of the restframe UV spectra between age and metallicity (cf. 
above) no clear indication on the galaxian metallicity can be derived, in contrast to the 
claim of KESR. Good fits to the available data can even be found with solar metallicity 
starburst templates. 

5) Depending on the star formation history and age one may or may not expect 
intrinsic Lya emission, i.e. an important H n region around the object. The apparent 
absence of observed Lya emission docs therefore not provide much insight. 

The theoretical templates can also be used to estimate the stellar mass involved in the 
starburst or the star formation rate when constant star formation is assumed. For this 
aim we assume a typical redshift of z = 6.6, and the magnification fi — 25 determined 
by KESR. For constant SF we obtain SFR ~ (0.9 - 1.1) M yr" 1 (for a Salpeter IMF 
from 1 to 100 M Q ). For the best fit ages of ~ 400-570 Myr the total mass of stars formed 
would then correspond to ~ (3.6 — 6.3) x 10 8 M Q . The mass estimated from best fit burst 
models (of ages ~ 6-20 Myr) is slightly smaller, M± ~ (0.3 — 1) x 10 8 M Q . If we assume 
a Salpeter IMF with Mi ow =0.1 M Q the mass and SFR estimates would be higher by a 
factor 2.55, and in good agreement with the values derived by KESR and Egami et al.. 
In all the above cases the total luminosity (unlensed) is typically Lboi ~ 2 x 10 10 L Q . 



4. z > 6 starbursts: with Herschel and ALMA, and now . . . 

Let us now assume that starburst galaxies with dust exist at z > 6 and briefly examine 
their observability with facilities such as Herschel and ALMA. To do so we must assume 
a typical galaxy spectrum including the dust emission. For simplicity we here adopt the 
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SED model by Melchior et al. (2001) based on PEGASE.2 stellar modeling, on the Desert 
et al. (1990) dust model, and including also synchrotron emission. Their predicted SED 
for a galaxy with an SFR and/or total luminosity quite similar to that estimated above 
for HCM6A is shown in Fig. 3. 

Moderate starburst - L, R = 1 .8 x 10" L„ - SFR = 32 M.yr" 1 



Z=0.1 0.5 1 235 10 20 30 




L ~ 1.8 x 1O U L 



placed at redshift z =0.1, 0.5, 1, 2, 3, 5, 10, 20, and 30. The thresholds of the JWST (here 
NGST), PACS and SPIRE onboard Herschel, and ALMA are also presented. Figure taken from 
Melchior et al. (2001) with kind permission 

Figure 3 shows the exquisite sensitivity of ALMA in the various bands allowing in 
principle an easy detection of such objects up to redshift ~ 10 or even higher! 

On the other hand, with the sensitivity of PACS and SPIRE on Herschel blank held 
observations of such an object are limited to smaller redshift (z < 1-4). However, already 
with a source magnification of [i ~ 3-10 or more the "template galaxy" shown in Fig. 
3 becomes observable with SPIRE at ~ 200-670 fim. In fact, such magnifications (and 
even higher ones) are not exceptional in the central parts of massive lensing clusters. 
E.g. in our near-IR search conducted in two ISAAC helds (~ 2.5x2.5 arcmin 2 ) of two 
lensing clusters, a fair number (~ 10-20) of z > 6-7 galaxy candidates with fj, > 5 are 
found (Pello et al. 2004b and these proceedings, Richard et al. 2004, in preparation). 
More than half of them have actually magnihcations [i> 10. Such simple estimates show 
already quite clearly the potential of strong gravitational lensing to extend the horizon 
of SPIRE/Hcrschel observations beyond redshift z>5\ 

Obviously a more rigorous feasibility study must also address the following issues: How 
frequent is dust present in high-z galaxies? and up to what redshift? We now have some 
indications for dust in one lenscd z = 6.56 galaxy (see Section 2) and of course in high-z 
quasars. But how general/frequent is this? How typical is the SED adopted above? The 
long wavelength emission due to dust depends on various parameters such as metallicity, 
the dust/gas ratio, geometry, the ISM pressure etc. Furthermore spatial resolution and 
source confusion are key issues which must be addressed and which should vary quite 
strongly between blank helds and cluster environments. Last, but not least, the held 
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of view of the various instruments is determinant for the efficiency with which high-z 
candidates can be found and studied. Several of these issues have already been partly 
addressed earlier (cf. the 2000 Herschel conference proceedings of Pilbratt et al. 2001, 
also Blain et al. 2002). 

It is evident that various ground-based and space bourne facilities and instruments 
will be used together to provide an optimal coverage in wavelength, spatial resolution 
and field size, and to obtain imaging as well as spectroscopy. Near-IR wide field imagers 
and near-IR multi-object spectrographs on 8-10m class telescopes and later with ELTs 
will undoubtably "team up" with the JWST, Herschel and ALMA to explore the first 
galaxies in the Universe and their evolution from the Dark Ages to Cosmic Reionisation. 
The wonderful power offered by gravitational lensing will continue to provide deeper or 
"enhanced" views of prime interest for the exploration of the early Universe. 
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